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Abstract 

fc-submodular functions, introduced by Huber and Kolmogorov, are functions defined 
on {0, 1 , 2 ,..., fc}" satisfying certain submodular-type inequalities, fc-submodular functions 
typically arise as relaxations of NP-hard problems, and the relaxations by fc-submodular 
functions play key roles in design of efficient, approximation, or fixed-parameter tractable 
algorithms. Motivated by this, we consider the following problem: Given a function / : 
{1,2,..., fc}"- —RU{-|-oo}, determine whether / can be extended to a /c-submodular function 
g : {0,1,2,..., /c}" —)• R U {-l-oo}, where g is called a fc-submodular relaxation of /, i.e., the 
restriction of g on {1, 2,..., fc}” is equal to /. 

We give a characterization, in terms of polymorphisms, of the functions which admit a 
fc-submodular relaxation, and also give a combinatorial 0((fc")^)-time algorithm to find a 
fc-submodular relaxation or establish that a fc-submodular relaxation does not exist. Our 
algorithm has interesting properties: (1) If the input function is integer valued, then our 
algorithm outputs a half-integral relaxation, and (2) if the input function is binary, then our 
algorithm outputs the unique optimal relaxation. We present applications of our algorithm 
to valued constraint satisfaction problems. 

Keywords: /c-submodular function, fc-submodular relaxation, valued constraint 

satisfaction problems 


1 Introduction 


A k-submodular function (Huber and Kolmogorov [B]) is a function / on {0,1,2,..., satisfying 
the following inequalities, 

f{x) + f{y)> f{xny) + f{xUy) (x, y € {0,1, 2,... , A:}”), (1) 


where binary operations fl, U are dehned by 


{x n y)i 


Xi if Xi = yt, 
0 if Xi 7^ yi, 


{x U y)i 


/ 

Xi if Xi = yi, 

0 if 0 A Xi A yi 7^ 0, 
yi if Xi = 0, 

Xi if yi = 0 

\ 
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Project. 

^Department of Mathematical Informatics, Graduate School of Information Science and Technology, University 
of Tokyo, Tokyo, 113-8656, Japan. 

Email:{hirai, yuni_iwamasa}@mist.i.u-tokyo.ac.jp 


1 




for X = {xi,X 2 , ■ ■ ■ ,Xn) and y = {yi,y 2 , ■ ■ ■ ,yn)- Observe that 1-submodular functions are 
submodular functions and 2-submodular functions are bisubmodular functions (see [2]). 

A:-submodular functions typically arise as relaxations of NP-hard problems, and the relax¬ 
ations by /c-submodular functions, /c-submodular relaxations, play key roles in the design of 
efficient, approximation, or fixed-parameter tractable (FPT) algorithms. For a function / on 
{1,2,..., k}^, a k-submodular relaxation [H [9] of / is a function g on (0,1,2,..., A:}" such that 
g is fe-submodular and the restriction of g to {1, 2,... , A;}"' is equal to /. Gridchyn and Kol¬ 
mogorov [3] showed that the Potts energy function, a generalization of the objective of multiway 
cut, has a natural A:-submodular relaxation, and that this relaxation is useful in computer vision 
applications. Iwata, Wahlstrom, and Yoshida [9] developed a general framework of FPT algo¬ 
rithms with introducing the concept of a discrete relaxation, where a A:-submodular relaxation is 
a primary and important example of discrete relaxations. Hirai [5j introduced a class of discrete 
convex functions that can be locally relaxed to A:-submodular functions, and designed efficient 
algorithms for some classes of multiflow and network design problems. 

In view of these appearances and applications of A:-submodular functions, it is quite natural 
and fundamental to consider the following problem: Given a function / on {1, 2,..., k}'^, deter¬ 
mine whether there exists a A:-submodular relaxation of /, and find a A;-submodular relaxation 
if it exists. 

The main results of this paper are a characterization of those functions which admit k- 
submodular relaxations, and a fast combinatorial algorithm to find a A:-submodular relaxation. 
Let [A;] := {1,2,..., A:} and [0, A:] := [A:] U {0}. In this paper, functions can take the infinite 
value -|-oo, where a < -|-oo and a -|- oo = -|-oo for o G M. The A;-submodular inequality ([I]) is 
interpreted in this way. (In the case of f{x) = -|-oo or f{y) = -|-oo, ([I|) trivially holds even if 
f{x n y) = - 1-00 and f{x Li y) = -|-oo.) Let M := M U {-|-oo}. For a function / : H"' —M, let 
dom / := {x G \ f{x) < -|-oo}. We show that the A:-submodular extendability is characterized 
by a certain operation on [A:]. Let us define a ternary operation 9 : [A:]^ [k] by 


9{a, b, c) 


a if a = b, 
c if a ^ b. 


The ternary operation 6 is extended to a ternary operation ([A:]”)^ —)■ [A:]” by {9{x,y, z))i = 
6{xi,yi, Zi). Note that 0 is a majority operation (in the sense of [E]), since 6{a, a, b) = 9{a, b, a) = 
9(b, a, a) = a. In universal algebra, 6 is known as the dual discriminator [T]; this fact was pointed 
out by A. Krokhin and the referees. 

Theorem 1. A function f : [A;]" —)■ M admits a k-submodular relaxation if and only if6{x, y, z) G 
dom / for all x,y,z € dom /. 

Therefore the class of A;-submodular extendable functions is defined by a polymorphism 9, and 
hence is closed under expressive power (see m)- Also the A:-submodular extendability depends 
only on the domain of /, i.e., the labelings attaining finite value. In particular, if dom / is the 
whole set [A:]"' (i.e., / : [A:]"' —>■ M), then / always has a A;-submodular relaxation; this fact has 
been noticed by Gridchyn and Kolmogorov [H p. 2325], but their proof is not correct.^ 

Based on Theorem [H we can obviously test for the existence of a A:-submodular relaxation in 
0{{k"')^) time by going through all labelings x,y,z € dom / and testing for the closure under 9. 
However this method cannot find a A:-submodular relaxation even if it exists. We will present a 

^They claimed that a fc-submodular relaxation g of arbitrary / is obtained by setting g{x) = f{x) for x £ [fc]" 
and g{x) = C for x £ [0, fc]"' \ [fc]", where C < minj.g[j.]n f{x). This is not trne. Indeed, consider / : [2]^ —>■ R such 
that /(I, 2) 1 and f[x) := 0 for other x £ [2]^. Let g : [0, 2]^ —>■ R be defined by g{l, 2) := 1 and g{x) 0 for 

other X £ [0, 2]^. Then g is not fc-submodular since 1 = g(0, 0) + g{l, 2) > g{l, 0) + g(0, 2) = 0. 
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combinatorial 0((/c”)^)-time algorithm to find a /c-submodular relaxation. Our algorithm reveals 
interesting and unexpected properties of the space of fe-submodular relaxations: the existence 
of a half-integral /c-submodular relaxation and the existence of a unique maximal /c-submodular 
relaxation in the case of n = 2. 

Theorem 2. There exists an O -time algorithm to determine whether a function f : 

[/c]” ^ M has a k-submodular relaxation, and to construct a k-submodular relaxation g if it 
exists, where g has the following properties: 

1. If f is integer valued, then g is half-integer-valued. 

2. If n = 2, then for every k-submodular relaxation g' of f it holds that 

g{x) > g'{x) {x G dom g). 

Namely g is the unique maximal k-submodular relaxation of f. 

In particular, our algorithm outputs a half-integral and optimal /c-submodular relaxation if 
n = 2. This solves, in the special case of binary /c-submodular relaxations, a question raised 
by [9]: Is there a way to decide the existence of discrete relaxations in general? 

The rest of this paper is organized as follows. In Section [21 we present applications of our 
algorithm to valued constraint satisfaction problems (VCSPs), where we utilize a recent remark¬ 
able result by Thapper and Zivny m that /c-submodular VCSPs can be solved in polynomial 
time. (The oracle tractability of /c-submodular function minimization is one of the prominent 
open problems in the literature; see laio As a consequence of properties 1 and 2 in Theorem 
121 our algorithm always constructs a half-integral /c-submodular relaxation for integer-valued 
VCSPs, and an “almost best” /c-submodular relaxation for binary VCSPs. We also present an 
application to the maximization problem, where we utilize a recent result by Iwata, Tanigawa, 
and Yoshida [8] on /c-submodular function maximization. In Section [3l we prove Theorem [1] and 
121 Our algorithm is based on the Fourier-Motzkin elimination scheme for linear inequalities. We 
show that the system of /c-submodular inequalities has a certain nice elimination ordering, and 
the Fourier-Motzkin elimination can be greedily carried out. 

2 Application 

Our algorithm is useful in VCSPs. Let us introduce VCSPs briefly; see m for detail. Let D 
be a finite set, called a domain. By a cost function on D we mean a function / : D” —>■ M for 
some natural number r = rf, called the arity of f. A set of cost functions is called a language 
on D. For a language C, a pair (/, a) of f € C and a : {1,2,... , r^} — {1, 2,... , n} is called 
a constraint on C. An instance of VCSP over language C, denoted by VCSP(£), is a triple 
I = (n, D, C) of the number n of variables, domain D, and a finite set C of constraints on C. The 
task of VCSP(T) is to find x = (xi ,..., Xn) G that minimizes 

fi{x') . 'y ^ f ...,Xu(j.. 

Let OPT(I) := minx fi{x). 

In the case where D = [0, fc] and C consists of fc-submodular functions, we call VCSP(T) a 
k-submodular VCSP. Thapper and Zivny m proved the polynomial solvability of /c-submodular 
VCSPs (see [TU] for the journal version). 

Theorem 3 (see mm)- k -submodular VCSPs can be solved in polynomial time. 

A k-submodular relaxation of an instance I = (n, [k],C) is an instance P = (n, [0, k],C') such 
that C is obtained by replacing each cost function in C with its /c-submodular relaxation. Notice 
that /// is a /c-submodular relaxation of fj. 
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A:-submodular autarky. From a minimizer of a A:-submodular relaxation, we obtain an au¬ 
tarky, a partial assignment of variables that keeps OPT, on the basis of the following property 
(called persistency). 

Theorem 4 (see lai). Let f be a function on [fe]” and g a k-submodular relaxation of f. For 
any minimizer y* G [0, k]^ of g, there exists a minimizer x* G [k]^ of f such that x* = y* for all 
i with y* 7 ^ 0. 

By our algorithm, for an instance I = (n, [k],C), we can construct a /s-submodular relaxation 
r, if it exists, in 0(|C|(A:'’)^) time, where r is the maximum arity of a function in C. This is a 
polynomial time algorithm in VCSPs. By Theorem [3l we obtain an optimal solution y* of I' in 
polynomial time. By Theorem HI in solving I, we can fix Xi to y* for all i with y* ^ 0. This 
contributes to reducing the size of VCSP. 

FPT algorithm. Iwata, Wahlstrom, and Yoshida [9] present an application of A:-submodular 
relaxation for FPT algorithms. Suppose that L consists of integer-valued cost functions. For an 
instance 1 = (n, [k],C) of VCSP(/1) and a A:-submodular relaxation I' = (n, [0, A:],C') of /, the 
scaling factor of I' is the smallest integer c such that c ■ g is integer valued for all g ^ C. Let 
d := OPT(/) — OPT(I'). Then we can solve an instance I in polynomial time, provided k^'^ is 
fixed. 

Theorem 5 (see^ [9l Lemma 1]). Let C be a language on [k] consisting of integer-valued cost 
functions. Let I be an instance of VCSP(T), and I' a k-submodular relaxation of I with the 
scaling factor c. Let d := OPT(/) — OPT(/'). We can solve the instance I by solving a k- 
submodular VCSP at most k^'^ times. 

Our algorithm constructs a fe-submodular relaxation I' with c = 1 or 2, though we do not 
say anything about the magnitude of d in general. In the binary case (r = 2) that includes many 
important VCSPs, our relaxation is an almost best fc-submodular relaxation in the following 
sense: Our relaxation has the smallest d among all fc-submodular relaxations. Indeed, for any 
/c-submodular relaxation L", it holds that fr{x) > fi"{x) by property 2 in Theorem [21 and hence 
d < OPT(/) —OPT(/"). Moreover our relaxation L' has the smallest cd, except for the case where 
there exists another fc-submodular relaxation I” with scaling factor 1 and OPT(I) — OPT(/") < 
cd. In such a case, the obtained cd is still a 2-approximation. 

In FPT applications, d* = OPT(/) is a more desirable parameter than d = OPT(/) — 
OPT(/'), since d* depends only on input I (see [9]). If OPT(I') > 0, then d < d*, and we 
can use d* as an FPT parameter. This is in the case where each function in I' is nonnegative 
valued. This leads to the notion of a nonnegative (-valued) k-submodular relaxation. For binary 
functions, our algorithm returns a nonnegative /c-submodular relaxation if it exists. This fact will 
be useful in design of FPT algorithms for binary VCSPs. It should be noted that /c-submodular 
relaxations for special binary functions (given in [21 Section 4.1]), are the same as relaxations 
obtained by our algorithm. 

Maximization. Nonnegative /c-submodular relaxation also has a potential to provide a unified 
approach to maximization. We here consider maximization of functions having no -|-oo. Following 
mm. just recently, Iwata, Tanigawa, and Yoshida [8] presented a 1/2-approximation algorithm 
for nonnegative /c-submodular function maximization. 

*Note that our definition of fc-submodular relaxation is slightly different from the one given by [9], where the 
definition in [9] requires one more condition miny = min/. Theorem [5] holds in our setting since the proof does 
not use this condition. 
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Theorem 6 (see [ 8 l Theorem 2.3]). There exists a polynomial time randomized 1/2-approximation 
algorithm for maximizing nonnegative k-submodular functions (given by value oracle). 

The maximum of a fe-submodular function g on [0, A:]” is always attained at [A:]"' (see [ 8 l 
Proposition 2.1]). Indeed, for any maximizer 2 G [0,A:]”, choose any zi,Z 2 G [A:]” with 2 = 
zi n Z 2 = zi U Z 2 \ both zi and Z 2 are maximizers by A:-submodularity, g{zi) + g{z 2 ) > 2,g{z). 
Therefore if we have a nonnegative A:-submodular relaxation of given / and its value oracle, then 
by using Iwata-Tanigawa-Yoshida algorithm we obtain an approximate maximum solution of 
/ with factor 1/2 on average. Our algorithm is again useful for the case where / is given as a 
VCSP form, i.e., the sum of small arity functions. If we obtain a nonnegative relaxation for each 
summand, then this approximation scheme is applicable. 

Further study on nonnegative A:-submodular relaxation is left to future work. 


3 Proofs 

For X G [0,A:]”, let Z{x) denote the number of indices i with Xi = 0. For A C \k]^, let Ce{A) 
denote the minimum subset X of [A;]” containing A such that 9{x,y,z) G X for all x,y,z G X. 
For B C [0, A:]”, let Cn{B) (resp., C'n,u(.B)) denote the minimum subset X of [0, k]"^ containing B 
such that xHy G X (resp., xHy, xUy G X) for all x,y & X. Note that all sets, C' 5 i(-), Cn(-)) C'n,u(‘)) 
are uniquely determined. In particular, A C 0 (A), B i—>• Cn{B) and B i-)- Cn,u{B) are closure 
operators. Observe that 8 can be represented by U as follows: 

0 ix,y,z) = {{xUy)Uz)U{xUy). (2) 

3.1 Proof of Theorem [T] 

Lemma 7. For all A C [A:]”, it holds that Cu^u{Ce{A)) = Cn(Ce(^)). 

Proof. The inclusion (D) is obvious. Therefore it suffices to prove that Cn(C'e(A)) is closed under 
U. Take arbitrary x,y G Cn{C0{A)). Our goal is to show xUy G Cn(C'e(A)). By the dehnition 
of Cn, there are ..., ..., y* G Ce{A) such that x = x^ FI • • • FI x^ and y = y^ FI • • • FI y* 

(note that FI is associative). Let X := {i j Xj / 0} and Y := {j \ yj 0}. First we show that 
there exists u G Ce^A) such that Ui = x, for i G X and Ui = yi foic i G Y \ X . By the dehnition 
of U, X, and Y, we have 


{xUy)i 


0 if Xj = y* = 0 or 0 / Xi / yi 7 ^ 0 , 
< Xi if i G X \ y or Xi = yi 7 ^ 0, 

,yi ifiGy\X. 


Let y \ X = {ji, • • • ,ja}- For all j G X \ X, there exists a pair {pj,qj) of indices in {1, 2,.., s} 
such that Xj-’ 7 ^ x^^, since Xj = 0. Let := 9{xPpx^Py^). Then we have u/ = x^^ = x^-’ = Xi 
for i G X, and = yj = yj. Dehne uB-'-F 9 {uB--dk-i {2 <k < a). It is easily seen 

that uB-B £ C 0 (A), = Xi for i G X, and = y* for i G y\X. Similarly, there exists 

V G Cg{A) such that Vi = yi lor i G Y and Xj = Xi for i G X \ Y. Hence uH v G C'n(C'e(H)). 
It holds that {u FI v)i = (x U y)i for all i G X U X. Therefore the set B := {z G Cn{Cg{A)) \ 
Zi = (x U y)i for alH G X U y} is nonempty. 

Take z G B with maximum Z{z). We show z = x U y (implying x U y G C'n(C'e(H)), 
as required). Suppose to the contrary that z 7^ x U y. By assumption, there exists I such that 
zi ^ xi = yi = 0. For this I, there exist p, q, r such that xf 7^ xf and y[ 7^ zp since xi = yi = 0. Let 

:= 9{xP,x'^,y^). It holds that w/ = Xi for i G X, and w/ = y[ 7 ^ zp Define w := 9{w^,u,y^) 
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{u G C 0 {A) such that Ui = Xi ior i € X and Ui = -jji ior i € Y \ X). It is clear that Wi = Xi for 
i € X, tCj = y* for i € y \ X, and wi = y[ 7 ^ zi. Therefore zPtw ^ B and Z{z) < Z{zr\ w), since 
{z n w)i = 0. However this is a contradiction to the maximality of z. Thus z = xUy. □ 

Lemma 8. For all A C [k]'^, it holds that Cn,u(^) H [k]"^ = Cg{A). 

Proof. By ([2]), we have C'n,u(^) 2 Co{A) D A. Since C'n,u is a closure operator, it holds that 
Cn,u{A) = Cn,u(C'n,u(^)) 2 C\^^u{Ce{A)) 3 C'n,u(^)- In particular, C'n,u(C' 6 »(^)) = C'n,u(^)- 
Hence C'n,u(>I)n[fc]" = C'n,u(C' 6 i(^))n[A:]"-. By LemmalU Cn,u{Ce{A))P[kY = Cn(C' 6 i(^))n[fc]”. 
Here Cn{Cg{A))P[kY' = Cg(A) holds. Indeed, for x G C'n(C'e(^))n[A:]”, there are x^,x‘^, ..., G 
Cg{A) with X = n • • • n x™’. Then x^ = x^ = • • • = x^ = x must hold, otherwise x includes 

O, contradicting x G [fc]”. Thus x G Co{A) and Cn{Cg{A)) n [/c]” = Ce{A). Consequently, 

Cn,u{A) n [kf = Ce{A). □ 

Proposition 9. A function f : [k]^ —>■ M admits a k-submodular relaxation if and only if 

(^□.□(dom /) n [k]'^ = dom /. 

Proof (only-if part). Suppose that / has a fc-submodular relaxation g. By the equation dom / = 
dom y n [A:]”, we have 

C'n,u(dom /) D [A:]"" ^ dom / = dom g D [A:]"" = C'n,u(dom g) D [A:]"" 

2 (Cn,u(dom g n [A:]'")) fl [A:]” = C'n,u(dom /) n [A:]”, 

where dom g = C'n,u(dom g) follows from the fact that the domain of a A:-submodular function 
is closed under □, U. Thus C'n,u(dom /) n [A:]” = dom /, as required. 

(if part). Assume C'n,u(dom /) n [A:]"’ = dom /. Let N := |C'n,u(dom /)| and N' := |dom f\. 
We can consider a function g : [0, A:]"’ —)■ M with dom g = C'n,u(dom /) as a vector g G , 
where the xth component Qx of g for x G dom g is defined as g{x). By definition, the set of 
all A:-submodular functions is defined by linear inequalities ([T]), and hence forms a polyhedron 
P = {y G I Ag < 0}. Therefore, the set of functions which admit a A:-submodular relaxation 
can be considered as the projection P' := {/ G | y G P, / is the projection of g to of 

P. Let us prove that the projection P' is equal to . 

We can obtain P' by the elimination of all variables gx (x G dom g \ dom /) by using the 
Fourier-Motzkin elimination method. We repeatedly eliminate variables gx by taking an index 
X with maximum Z[x), i.e., the number of zeros in x, in each step. Suppose that an index x is 
chosen in the first step. Then coefficients of gx in Ag < 0 are positive or zero. Indeed, assume 
that there exists an inequality such that the coefficient of gx is negative. Namely, there is a 
nontrivial inequality 

g{x n y) + y(x U y) - y(x) - y(y) < 0 (3) 

with xriyT^xy^xLly. By the definition of x, we have Z[x) > y(x □ y). By the definition of 
n, we have Z{x) < Z{x □ y). Thus Z{x) = Z{x □ y), x = x □ y, and y = x U y. This means 
that the coefficient of gx is positive in all inequalities containing gx- So the linear inequality 
system of the projection is obtained by simply removing all inequalities containing variable gx- 
Now suppose that the set S of variables has been eliminated by the Fourier-Motzkin method, 
and the corresponding system of inequalities consists of the original inequalities not containing 
variables in S, as above. In the next step, the Fourier-Motzkin procedure chooses an index x 
with maximum Z(x) over dom g\S. By a similar argument, the coefficients of gx in the current 
inequalities are positive or zero; inequalities having negative coefficients at gx have already been 
removed in the previous steps. Thus there are finally no inequalities. (Note that there exist 
no inequalities only containing elements in dom /, since if x, y G dom / and x 7 ^ y, then 
xny = xUyG dom g \ dom /.) This means that P' = , as required. □ 
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By the proof of Proposition [9l the A:-submodular inequalities can be represented as follows: 
For any G dom g such that z = xFly and Z{z) > max{Z(x), Z(y)}, 

^(^) < ( 5 ( 2 ;)+ 5 ( 2 /)) ifxny = xUy, 

I g{x) + g{y) — g{x U y) otherwise. 

We are now ready to prove Theorem[TJ By Lemma[8]and Proposition O a function / : [A:]” —>■ 
M admits a fc-submodular relaxation if and only if (^^(dom /) = dom /. It is clear that this 
statement is the same as Theorem [TJ 

3.2 Proof of Theorem [2] 

We present an algorithm with the claimed properties in Algorithm [TJ Let us prove that Algo- 


Algorithm 1 A:-submodular relaxation 
Initialize g^^'^ as follows: 




if X G dom /, 
otherwise. 


for i = 1 to n do 

■= 

for all x,y £ dom such that Z{xny) = i do 

if x n y = x U y then 

if 0 ^*) (x n y) > (x) -|- (y)) 12 then 

y«(xny) := (x) + (y))/2 

end if 

else if X U y 0 dom then 

retnrn / has no /c-submodular relaxation, 
else 

if y^*^ (x n y) > (x) -|- (y) — (x U y) then 

y W (x n y) := (x) -|- (y) — (x U y) 

end if 
end if 
end for 
end for 
return y^”^ 


rithm [T] correctly determines whether a function / : [A:]"' —?> M has a A:-submodular relaxation, 
and constructs a A:-submodular relaxation if it exists. First we show that if Algorithm [T] returns 
“/ has no A:-submodular relaxation,” then the input function / actually has no A:-submodular 
relaxation. To prove this statement, we show the following claim. 

Claim 10. It holds that dom y^*) = {z £ Cn(dom /) | Z{z) <!}. 

Proof of Claim[TR The inclusion (C) is obvious. We prove (A) by induction on i. The case 
A = 1 is trivial. Indeed, any element 2; of dom y = C'n(dom /) with Z{z) = 1 can be written 
as z = X n y with x, y G dom /. For all 2; G C'n(dom /) \ dom / with Z{z) = i -|- 1, there exist 
x,y G C'n(dom /) such that max{Z(x), Z{y)} < i and 2; = x PI y. By the induction hypothesis, 
we have x,y G dom g^^\ and hence 2 G dom This completes the induction step. □ 
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Here we consider the case of returning “/ has no /c-submodular relaxation.” In this case, 
for some step i > 1, there are € dom such that Z{x □ y) = i, x n y ^ x U and 

X U y 0 dom Then x,y € dom and Z{xr\y) > Z{xUy). Therefore xUy ^ 

dom = {z G C'n(dom /) | Z{z) < i — 1} (by Claim fT(l . On the other hand, it is obvious 

that xUy G {z G C'n,u(dom /) | Z{z) < i — 1}. This means that C'n(dom /) / C'n,u(dom /). By 
Lemma[71 it necessarily holds that C 6 i(dom /) 7 ^ dom /. By Theorem[Tl there is no fc-submodular 
relaxation for /. 

Next we consider the case of returning . Returning means Cn,u(dom /) = C'n(dom /) 
since for all x, y G C'n(dom /), x U y G Cn(dom /). Therefore Cn,u(dom /) 0 [/c]"' = dom /, and 
/ admits a fc-submodular relaxation by Proposition [9j Furthermore defined by Algorithm 1 
satisfies ([H). Thus is a /c-submodular function. 

Next we show the half-integrality property (property 1 in Theorem [2]). 

Proposition 11. Suppose that a function f : [k]"' —>■ M has a k-submodular relaxation. Let g he 
a k-submodular relaxation of f constructed by Algorithmic For all z G dom g \ dom f, there 
exist X, y G dom g with z = xHy and Z{z) > max{Z(x), Z{y)} satisfying 1 or 2: 

1. x,y G dom / and g{z) = ^(y(x) +g{y)); 

2. xUyi^xUy and g{z) = g{x) -F y(y) - y(x U y). 

In particular, if f is integer valued, then g is half-integer-valued 

Proof. We will prove this by induction on Z{z). By Algorithm [H for some x,y G dom g with 
z = xHy and Z{z) > max{Z(x), Z{y)}, the value g{z) is equal to (y(x) -\-g{y))/2 if xFly = xUy, 
and y(x) -|- y(y) — g{x U y) otherwise. If Z{z) = 1, then Z{x) = Z(y) = 0 implying x, y G dom /, 
X n y = X U y, and g{z) = (y(x) -|- y(y))/2; we are in 1. Suppose that Z{z) > 2, and that x,y 
satisfy neither 1 nor 2. Then x □ y = x U y, g{z) = (y(x) -|- y(y))/2, and therefore one of x, y does 
not belong to dom /. Say x 0 dom /. Here 1 < Z{x) < Z{z) = Z{x n y). So by the induction 
hypothesis applied to x, we only consider two cases: 

Case 1: y(x) = -(y(x^) -|- y(x^)) for some x^,x^ G dom / with x^ FI x^ = x; 

Case 2: y(x) = y(x^) -|- y(x^) — y(x^ U x^) for some x^, x^ G dom g with x^ U x^ 7 ^ x^ □ x^ = x 
and Z{x) > max{Z(x^), Z(x^)}. 

Let Z := {i \ Zi ^ 0}. Let X := {i ^ Z \ Xi ^ 0} and y := {i ^ Z \ yi ^ 0}. Note that 
A = T, 0 7 ^ Xj 7 ^ yi 7 ^ 0 for / G A, and Xj = y^ 7 ^ 0 for z G Z, since z = xny = xLly. Hence 
Z{x) = Z{y) > 1, i.e., y G dom g \ dom /. 

In Case 1, it holds that Xj = xj = x? for i G A U Z and x^ 7 ^ x? for i 0 A U Z. We obtain 


9 {z) = ]^{g{x)+g{y)) (5) 

= \{^9{x^) + \9{x^'^ + g{y'^ ( 6 ) 

> ^ (^\9{x^ n y) -h ^y(x^ U y) -h ^g{x‘^ n y) -b ^g{x‘^ U y)^ ( 7 ) 

= \ ( 5 ( 2 ) + ]^9{x^ U y) + iy(x2 U y)^ (8) 

>]^{ 9 {z)+ 9 {z))=g{z). (9) 
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Indeed, ([5]) = ([6]) follows from the assumption of Case 1, and (l6|) > ([7|) follows from the k- 
submodularity {g{x^) + g{y) > g{x^ n y) + g{x^ U y) and + g{y) > g{x‘^ n y) + y(x^ U y)). 

Since Xi = xj = xf = yi for i £ Z, Xi = xj = xf ^ yi for i £ y, and yi = 0 for i ^ y U Z, 
it holds that n y = n y = z. Hence dZ]) = ([8]). Since (x^ U y)i = (x^ U y)i for i £ Z, 
(x^Uy)j = (x^Uy)j = 0 for i € Al, and x\ = (x^Lly)i ^ (x^Uy)j = x‘f foi i ^ XUZ, it holds that 
(x^ U y) n (x^ Uy) = (x^ U y) U (x^ Uy) = z. Hence ([8|) > ([9]) follows from the /c-submodularity. 
This means that all inequalities are equalities. Therefore y(x^) + y(y) = y(x^ n y) + y(x^ U y) by 
dSD = (0. (It is also true that y(x^)+y(y) = y(x^ n y)+y(x^ U y).) Here y(x^ □ y) = y( 2 ;). Thus 
g{z) = y(x^) + y(y) — y(x^ U y). Moreover it holds that z ^ x^ \Jy. Indeed, for i 0 T U Z we 
have Zi = (x^ □ y)j = 0 and (x^ U y)j = xj ^ 0, since y^ = 0. So g{z) = y(x^) + y(y) — g{x^ U y) 
means that x^,y satisfy 2. 

In Case 2, it holds that Xj = x^ = x? for z € T” U ^ and (x^ □ x^)j = 0 for z 0 T” U We 


obtain 

9{z) = ^{9{x) + g{y)) (10) 

= ^(.9{x^) + 9{x^) - 9{x^ U x^) + y(y)) (11) 

> \{9{x^ n y) + y(x^ U y) + y(x^) - y(x^ U x^)) (12) 

> \i9{z) + 9{{x^ U y) n x^) + y((x^ U y) U x^) - y(x^ U x^)) (13) 

= \i9{z) + 9{z) + 9{x^ U x^) - y(x^ U x^)) = g{z). (14) 


Indeed, (fTU|l = (fTTT) follows from the assumption of Case 2, and (fTT]l > (fT^ follows from the 
fc-submodularity. Since Xj = xj = yi for i £ Z, Xi = x\ ^ yi for i £ X, and y* = 0 for z ^ T U .E, 
it holds that x^Hy = z. Hence (1121) > m follows from the /c-submodularity. Since (x^Uy)^ = xf 
for i £ Z, (x^ U y)j = 0 for i £ X, and (x^ U y)j = x) for z 0 T U E, it holds that (x^ U y) □ x^ = z 
and (x^ U y) U x^ = x^ U x^. Hence (jl3p = (jl4l) . This means that all inequalities are equalities. 
Therefore y(x^) + y(y) = g{x^ n y) + y(x^ U y) by (fTTI) = ([12]). Here y(x^ n y) = g{z). Thus 
g{z) = y(x^) + y(y) — y(x^ Uy)- Moreover it holds that z x^Uy. Indeed, there exists z ^ TUE 
such that xl ^ 0, since E(y) = E(x) > Z{x^). Hence Zi = (x^ □ y)j = 0 and (x^ U y)j = x\ ^ 0. 
So g{z) = y (x^) + 5 (y) “ 9 {x^ U y) means that x^, y satisfy 2. 

Next let us prove that if / is integer valued, then for all z £ dom g \ dom /, g{z) is half¬ 
integral. The proof is by induction on Z{z). The case Z{z) = 1 is trivial. Suppose that Z{z) > 2. 
For all z £ dom g \ dom /, there exist x,y £ dom g with z = x FI y and Z{z) > max{E(x), E(y)} 
satisfying 1 or 2: 

1 

1. x,y £ dom / and g{z) = -(y(x) +y(y)); 

2. X n y / X U y and g{z) = g{x) + g{y) - g{x U y). 

In case 1, it holds that g{z) is half-integral, since both y(x) and y(y) are integral. In case 2, 
it also holds that g{z) is half-integral, since y(x), y(y), and y(x U y) are half-integral by the 
induction hypothesis. Hence g{z) is half-integral, as required. □ 

Finally we establish property 2 in Theorem |21 

Proposition 12. Suppose that f : [k]"^ —>■ M has a k-submodular relaxation. Let g be a k- 
submodular relaxation of f constructed by Algorithm\J[ For every k-submodular relaxation g' of 
f, it holds that 

g{z) > g'(z) (z £ dom g). 
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Proof. Let be a set of /c-submodular relaxations of /. By Algorithm [U for any 2; € dom g \ 
dom /, there exist three cases of x, y G [0, k]'^ with z = xUy m. defining g{z) as follows: 

Case 1: Z{x) = Z{y) = 0, 

Case 2: Z{x) = Z{y) = 1, 

Case 3: Z{x) = l,Z{y) = 0. 

By Proposition [m we do not have to consider the case of g{z) = {g{x) + g{y))/2 with x G 
dom g \ dom /. So we can represent the three cases as follows: 

Case 1': g{z) = ^{g{xi,X2) + 5(2/1,2/2)) for some xi,X2,2/1,2/2 G [k], 

Case 2 ': g{z) = g{0, 0 ) = g{x, 0 ) + 5(0, y) - g{x, y) for some x, 5 G [k], 

Case 3 ': g{z) = 5(0,0) = 5 (a^i, 0 ) + 5(51,52) -5(0,52) for some xi,yi,y 2 e [k] with xi 7^ 51. 

It suffices to prove that for all g' G G, we have g(z) > g'{z) for all cases. Note that for 
X G dom / = dom g fl [k]'^ = dom g' fl it holds that g{x) = g'{x). Furthermore it holds that 
dom g' D C'n,u(dom /) = dom g, since dom g' D dom / and g' is A;-submodular. 

First we shall consider Case 1'. For any g' G G, it holds that g'(z) < (g'(xi,X 2 )-l-g'(yi,y 2 ))/^ = 
(g(xi,X2) + g(yi,y 2))/2 = g(z). So g'(z) < g(z), as reqnired. 

Next we consider Case 2 '. For any g' G G, it holds that 5'( 0 , 0 ) < 5'(x,0) +5^(0,5) — g'{x,y). 
In addition, we have 5'(x,0) < 5(x,0) and g'{0,y) < 5(0,5), since 5(x,0) and 5(0,5) are in 
Case 1 '. Hence 5'(0,0) < g'{x,0) +5'(0,5) - g'{x,y) < 5(^,0) +5(0,5) - 9 {x,y) = 5(0,0). So 
5^(0,0) < 5(0,0), as required. 

Finally we consider Case 3 '. We show that this case can reduce to Case 2' . It holds that 
5(xi,0) + 5(51,52) -5(0,52) > 5(3^1, 0 ) +5(0,52) - 5 (xi, 52 ) since 

(5(xi,0) +5(51,52) - 5(0,52)) - {g{xi,0) +5(0,52) - 5(xi,52)) 

= 5(51,52) + 5(3^1,52) - 25(0,52) > 0. 


Hence we have 


5(0,0) = 5(aii, 0 ) + 5(51,52) - 5(0,52) 

> 5(2^1,0) +5(0,52) - 5 (xi, 52 ) 

> 5 ( 0 , 0 ). 

So all inequalities are equalities. This means that 5(0,0) = 5(xi,0) + 5(0,52) — g{xi,y 2 ), and 
hence Case 3' can rednce to Case 2 ', as required. □ 

We found, by computer experiments, functions on [A:]” for k,n > 3 or k = 2,n > 4 such that 
the maximal relaxation does not exist, and the proposed algorithm does not ontput an optimal 
relaxation 5 , i.e., the minimum value of 5 is greater than the minimum value of any relaxation. 
(We do not know the case for k = 2,n = 3.) 
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